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We report on a search for gravitational waves from coalescing compact binaries, of total mass between 2 
and 35 Mq, using LIGO observations between November 14, 2006 and May 18, 2007. No gravitational-wave 
signals were detected. We report upper limits on the rate of compact binary coalescence as a function of total 
mass. The LIGO cumulative 90%-confidence rate upper limits of the binary coalescence of neutron stars, black 
holes and black hole-neutron star systems are 1.4 x 10" 2 , 7.3 x 10" 4 and 3.6 x 10~ 3 yr 1 L 10 respectively, 
where Lio is 10 10 times the blue solar luminosity. 

PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.60.Jd, 97.60.Lf, 97.80. -d 



In November 2005 the three first-generation detectors of the 
Laser Interferometer Gravitational-wave Observatory (LIGO) 
reached design sensitivity and began a two-year period of 
observations (known as the fifth science run, or S5) which 
concluded in October 2007 |Q]]. One of the most promising 
sources of gravitational-waves for LIGO is a compact binary 
coalescence (CBC); the inspiral and merger of binary neu- 
tron stars (BNS), binary black holes (BBH), or a black hole- 
neutron star binary (BHNS) H H H S H|j. These systems 
spiral together as they emit energy in the form of gravita- 
tional waves, finally merging to form a single object, which 
then settles down to equilibrium. Ground-based gravitational- 
wave detectors are most sensitive to waves with frequencies 
between ~ 40 and 1000 Hz, corresponding to the late stages 
of inspiral and merger. In this paper we report the results of 
search for gravitational-waves from binaries with total mass 
between 2 and 35 Mq and a minimum component mass of 
1 Mq in LIGO observations between November 14, 2006 and 
May 18, 2007. The results of a search for these systems in 
data taken from November 4, 2005 to November 14, 2006 
were reported in Ref. [7]. From May-October 2007, the Virgo 
gravitational-wave detector operated in coincidence with the 
LIGO detectors |Ht] and the LIGO data from that period are 
being analyzed together with the Virgo data. The joint analy- 
sis requires significant modifications to our analysis pipeline: 
therefore results of that search will be reported in a subse- 
quent publication. In contrast, the results presented here were 
obtained with substantially the same analysis pipeline used in 
Ref. 0| . 

No gravitational-wave signals were observed during this 
search and so we report upper limits on CBC rates using the 
upper limits of Ref. H as prior rate distributions. We summa- 
rize the analysis procedure and we present the search results 
and upper limits on CBC rates derived from LIGO observa- 
tions in the period November 4, 2005 to May 18, 2007. 



ThedData analysis pipeline: The data-analysis pipeline 
used in this search is fundamentally the same as that of 
Ref. 17[], thus here we only describe the major components 
and highlight differences to the previous search, referring to 
Refs. (a, ED for details. The most substantial change in this 
analysis is a modification to the way in which the significance 
of candidate events is compared to instrumental noise back- 
ground. In previous searches, the noise background was com- 
puted using the entire observation period by introducing an 
artificial time shift between data recorded at the two LIGO ob- 
servatories. The observation period is split into six four-week 
segments and one 18 day segment (referred to as "months") 
and the instrumental background is measured independently 
in each month, as the detector behavior varied over the course 
of the S5 run. Candidate triggers are therefore compared to 
a background that better reflects the instrumental behavior at 
the time of the candidate. Each month was searched indepen- 
dently for gravitational-wave candidates and in the absence 
of detections, the results from the months are combined (to- 
gether with the results from Ref. 1211) to set an upper limit on 
the CBC rate. 

We search for gravitational-wave signals when at least two 
of the LIGO detectors were operational. This comprised a 
total of 0.28 yr when all three detectors (the 4 and 2 km Han- 
ford detectors, denoted HI and H2, respectively, and the 4 km 
Livingston detector, denoted LI) were operational (H1H2L1 
coincident data), 0.10 yr of H1H2 coincident data, 0.02 yr of 
H1L1 coincident data, and 0.01 yr of H2L1 coincident data. 
Noise correlations between the colocated HI and H2 detectors 
cause our method of estimating the instrumental background 
using time-shifted data to fail, and so we do not search data 
when only the H1H2 detectors are operating. Approximately 
10% of data is designated playground and used for tuning our 
search pipeline. 

Post-Newtonian (PN) theory provides accurate models of 
the inspiral waveform predicted by general relativity up to the 
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innermost stable circular orbit (ISCO) H [M [H [H [H [3 
[l6ll . The frequency of the waveform from low mass bi- 
naries targeted in this search sweeps across the sensitive band 
of the LIGO detectors. Therefore, we search for signals from 
our target sources by match filtering the data with PN tem- 
plates terminated at ISCO. This method is suboptimal if a 
true signal differs from our template family due to unforeseen 
physical effects. Matter effects in BNS and BHNS are not 
included in our templates, but are expected to be important 
only at hig her frequencies lfl7l [Till . We construct template 
banks IU9I1 of restricted second order PN waveforms in the 
frequency domain fiol l20l l2lll such that no more than 3% of 
the signal-to-noise ratio (SNR) is lost due to the discreteness 
of the bank [22]. A "trigger" is generated if the matched-filter 
SNR of the strain data filtered against the template exceeds a 
threshold of 5.5 l23ll . We demand that triggers are coincident 
in time of arrival and mass l24ll in at least two of the three 
LIGO detectors. When all three detector are operating we 
can obtain (in principle) four possible types of coincidence: 
H1H2L1 triple coincident triggers and three different double 
coincident types: H1H2, H1L1 and H2L1. We discard H1H2 
double coincident triggers, due to the problems estimating the 
background for these triggers and discard H2L1 triggers when 
the HI detector is operating nominally (since the 4 km HI de- 
tector is more sensitive than the 2 km H2 detector). Coinci- 
dent triggers are subjected to consistency checks using signal- 
based vetoes ll25l 12a, l27ll . Times of poor detector data qual- 
ity are flagged using environmental and auxiliary data; trig- 
gers from these times are also vetoed 0. We construct two 
categories of data-quality vetoes depending on the severity of 
the instrumental artifact being flagged. In our primary search 
and upper limit computation we veto coincident triggers that 
fall in times from either category. We also consider detec- 
tion candidates in data with only the most severe category ap- 
plied in case a loud signal is present that may otherwise be 
vetoed. Surviving triggers are clustered in time and ranked by 
an effective SNR statistic, which is computed from the trig- 
ger's matched-filter SNR and the value of the \ 2 signal-based 
veto for that trigger J(J. After discarding playground data and 
times in both veto categories, a total of 0.21 yr of triple coin- 
cident data (H1H2L1) ,0.02 yr of H1L1 coincident data, and 
0.01 yr of H2L1 coincident data remain. In the absence of a 
detection, these data are used to compute upper limits on the 
CBC rate. 

The rate of instrumental noise artifacts is measured by ime- 
shifting data from the Livingston and Hanford observatories 
(HI and H2 data are kept fixed with respect to each other). 
The data are offset by more than the light-travel time between 
observatories, thus triggers which survive the pipeline are due 
to noise alone. We performed 100 such time-shifts to obtain 
a good estimate of the noise background in our search. CBC 
signals of higher-mass contain fewer gravitational-wave cy- 
cles in the sensitive band of our detectors: our signal-based 
vetoes are not as powerful. High-mass templates are therefore 
more sensitive to nonstationary noise transients and hence our 
false alarm rate (FAR) for these system is larger. In order 
to account for this mass-dependent behavior we compute the 
background for three different mass regions and compare fore- 
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FIG. 1: The posterior distribution for the rate of BNS coalescences. 
The dashed black curve shows the rate computed in Ref. @]. The 
solid black curve shows the result of this search using the previous 
analysis as a prior. The figure also shows the rate distributions for 
two of the individual months computed using a uniform prior. The 
improvement from month to month 5 is due to increasing detector 
sensitivity during this search. 

ground and background within each of these ranges. Specifi- 
cally, in each region we count the number of background trig- 
gers with effective SNR greater than or equal to a given fore- 
ground trigger; dividing this number by the amount of back- 
ground time analyzed gives us the FAR for that trigger. This 
allows us to define a single detection statistic for every trigger 
in each of the mass categories. The FAR can then be directly 
compared to obtain a ranking of the significance of the trig- 
gers, regardless of their mass 0]. 

Search results: The seven months of data were ana- 
lyzed separately using the procedure described above. No 
gravitational-wave candidates were observed with a FAR sig- 
nificantly above those expected from the noise background. 
The loudest trigger in this search was a triple coincident trig- 
ger with a FAR of 6 per year. This is consistent with the ex- 
pected background, since we searched 0.21 yr of data. The 
second and third loudest triggers had FAR values of 10 and 11 
per year respectively. Although we did not have any detection 
candidates, we exercised our follow-up procedures by exam- 
ining any triggers with a FAR of less than 50 per year. This 
exercise prepares us for future detections and often identifies 
areas where our search pipeline can be improved to exclude 
noise transients. 

In the absence of detection candidates, we use our observa- 
tions to set an upper hmit on the CBC rate. We follow the 
procedure described in IT281 129I [30ll and use the results re- 
ported in Ref. 0] as prior information on the rates. We present 
five different classes of upper limits. The first three limits 
are placed on binaries of neutron stars and/or black holes as- 
suming canonical mass distributions for BNS [mi = m.2 = 
(1.35 ± 0.04) M ], BBH [mi = m 2 = (5 ± 1) M Q ], and 
BHNS [mi = (5 ± 1) M Q , m 2 = (1.35 ± 0.04) M ] sys- 
tems. We also present upper limits as a function of the total 
mass of the binary and, for BHNS binaries, as a function of 
the black hole mass. We combine the results from each of the 
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seven months, along with the prior results from the first year 
analysis, in a Bayesian manner, using the same procedure as 
described in J3]. 

We first calculate upper limits on BNS, BBH and BHNS 
systems assuming the objects have no spin, and summarize 
the results Tables HI and iHl The rate of binary coalescences in 
a galaxy is expected to be proportional to the blue light lumi- 
nosity of the galaxy Oil] . Therefore, we place limits on the 
rate per Lio per year, where Lio is 10 10 times the blue solar 
luminosity (the Milky Way contains ~ 1.7Lio l32tl ). To cal- 
culate the search sensitivity, the analysis was repeated numer- 
ous times adding simulated signals with a range of masses, 
distance and other astrophysical parameters to the data. Ta- 
ble |n] shows the sensitivity of the LIGO detectors to coalesc- 
ing binaries quoted in terms of the horizon distance i.e., the 
distance at which an optimally oriented and located binary 
would produce an SNR of 8. There are a number of uncertain- 
ties which affect the upper limit calculation, including Monte 
Carlo statistics, detector calibration, distances and luminosi- 
ties of galaxies listed in the galaxy catalog {3l"1l and differ- 
ences between the PN templates used to evaluate efficiency 
of the search and the actual waveforms. The effect of these 
errors on the cumulative luminosity are summarized for the 
BNS search in Table U We marginalize over all of the uncer- 
tainties IT28ll to obtain a posterior distribution on the rate of 
binary coalescences. 

In Fig. [U we show the derived distribution of the rate of 
BNS coalescences. The distribution is peaked at zero rate 
because there are no detection candidates. We include the 
distribution for all searches previous to this one (which is 
our prior). In addition, we present the result that would be 
obtained from each month, were it analyzed independently 
of the others and of the previous searches. This provides 
an illustration of the amount that each month contributes to 
the final upper limit result and demonstrates the improve- 
ment in sensitivity of the detectors during the search. The 
upper limit is finally obtained by integrating the distribution 
from zero to TZ^o% so that 90% of the probability is con- 
tained in the interval. The results obtained in this way are 
^•90%,bns = 1-4 x lO^yr^LicT 1 ,^ 90 %,bbh = 7.3 x 
lO-V- 1 ^- 1 , and ^ 90%! bhns = 3.6xl0- 3 yr-iLio" 1 . 

Additionally we calculate the upper limit for BBH sys- 
tems as a function of the total mass of the binary, assum- 
ing a uniform distribution of the component masses. For 
BHNS systems, we construct an upper limit as a function of 
the black hole mass, assuming a fixed neutron star mass of 
«^ns = 1.35Af Q . These upper limits are shown in Fig[2] 

Finally, we present upper limits on coalescence rates where 
the spin of the components of the binary is taken into ac- 
count. Astrophysical observations of neutron stars indicate 
that their spins will not be large enough to have a signif- 
icant effect on the BNS waveform observed in the LIGO 
band 11331 13411 . Theoretical considerations limit the magni- 
tude of the spin, S, of a black hole to lie within the range 
< S < Gm 2 /c. However, the astrophysical distribution of 
black hole spins, and spin orientations, is not well constrained. 
Therefore, we provide a sample upper limit for spinning sys- 
tems using a spin magnitude and orientation distributed uni- 



Coincidence time 


H1H2L1 


H1L1 


H2L1 


Observation time (yr) 


0.21 


0.02 


0.01 


Cumulative luminosity (£10) 


490 


410 


110 


Calibration error 


23% 


23% 


26%) 


Monte Carlo error 


3% 


7% 


10% 


Waveform error 


31% 


32% 


31% 


Galaxy distance error 


16% 


16% 


3% 


Galaxy magnitude error 


19% 


19% 


17% 



TABLE I: Detailed results from the BNS search. The observation 
time is the time used in the upper limit analysis. The cumulative 
luminosity is the luminosity to which the search is sensitive above the 
loudest event for each coincidence time. The errors in this table are 
listed as one-sigma logarithmic error bars (expressed as percentages) 
in luminosity associated with each source error. 



Component masses 


1.35/1.35 


5.0/5.0 


5.0/1.35 


Aiorizon (Mpc) 


~ 30 


~ 100 


~ 60 


Cumulative lminosity 

(Lio) 


490 


11000 


2100 


Nonspinning upper 
limit (yr-^-Q 1 ) 


1.4 x 10~ 2 


7.3 x 10~ 4 


3.6 x 10~ :i 


Spinning upper limit 




9.0 x 10 -4 


4.4 x 10~ :i 



TABLE II: Overview of results from BNS, BBH and BHNS searches, 
^horizon is the horizon distance averaged over the time of the search. 
The cumulative luminosity is the luminosity to which the search is 
sensitive above the loudest event for times when all three LIGO de- 
tectors were operational. The first set of upper limits are those ob- 
tained for binaries with nonspinning components. The second set of 
upper limits are produced using black holes with a spin uniformly 
distributed between zero and the maximal value of Gm 2 jc. 

formly within the allowed values. This gives upper limits on 
the rate of BBH and BHNS systems of lZ 90 % yBBH = 9.0 x 

10~ 4 yr-^io" 1 andft 90%! BHNS = 4.4 xlO" 3 yr^LnT 1 . 
These rates are about 20% larger than the nonspinning rates. 

Discussion: We have searched for gravitational waves 
from CBCs with total mass between 2 and 35 Mq in LIGO 
observations between November 14, 2006 and May 18, 2007. 
No detection candidates with significance above that expected 
due to the background were found in the search. By combin- 
ing this search with our previous results, we set a new upper 
limit on the CBC rate in the local universe which is approxi- 
mately a factor of 3 lower than that reported in Ref. |7|]. This 
improvement is significant, even though we searched only two 
thirds as much data as in Ref. [7|]. It is due, in part, to improve- 
ments in detector sensitivity during S5 which increased the 
horizon distance. Moreover, the shorter analysis time and im- 
proved stationarity of the data, led to many of the months hav- 
ing a less significant loudest event than in the previous search. 
Both of these effects increased the luminosity to which the 
search was sensitive, thereby improving the upper limit. 

Astrophysical estimates for CBC rates depend on a num- 
ber of assumptions and unknown model parameters, and are 
still uncertain at present. In the simplest models, the coa- 
lescence rates should be proportional to the stellar birth rate 
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FIG. 2: The marginalized 90% rate upper limits as a function of 
mass. The upper plot shows limits for BBH systems as a function of 
the total mass of the system. The lower plot shows limits for BHNS 
systems as a function of the black hole mass, assuming a fixed neu- 
tron star mass of 1.35Mq . Here the upper limits are calculated using 
only H1H2L1 data since the relatively small amount of H1L1 and 
H2L1 data makes it difficult to evaluate the cumulative luminosity in 
the individual mass bins. 



in nearby spiral galaxies, which can be estimated from their 
blue luminosity 13 ill . The optimistic, upper end of t he p lau- 
sible rate range for BNS is 5 x 1CT 4 vr^I^Jpl |3rj|l and 
6 x 1(T 5 yr-^^ 1 for BBH and BHNS Q l38l The up- 
per limits reported here are ~ 1-2 orders of magnitude above 



the optimistic expected rates. With the next run starting in 
mid 2009, the Enhanced LIGO and Virgo detectors will be- 
gin operations with a factor of ~ 2 increase in horizon dis- 
tance. The total luminosity searched will increase by a factor 
of ~ 10, thereby bringing us close to the optimistic rates. The 
most confident BNS rate predictions are based on extrapola- 
tions from observed binary pulsars in our Galaxy' these yield 
realistic BNS rates of 5 x 10" 5 yr" 1 ^ 1 JHUBfl. Rate esti- 
mates for BBH and BHNS are less well constrained, but re- 
alistic estimates are 2 x 10 ~ 6 v r" 1 !^ 1 for BHNS $37^ and 
4 x 10~ 7 yr^Ljo 1 for BBH H. Thus, the expected rates are 
~ 2-3 orders of magnitude lower than the limits presented 
in this paper. The Advanced LIGO and Virgo detectors, cur- 
rently under construction, will increase our horizon distance 
by an order of magnitude or more, allowing us to measure the 
rate of CBCs in the Universe. 
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